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Abstract. The influence of the tip-substrate bias induced electric field in a scanning- 
tunneling-spectroscopy experiment on charged InAs nanocrystals is studied. Calculating the 
ground and first excited many-particle state for five electrons occupying the quantum dot 
reveals a Stark-induced reordering of states by increasing the electric field strength. It is shown 
that this reordering of states is accompanied by a symmetry change of the local density of states 
(LDOS), which in principal is observable in a wave-function mapping experiment. Since in the 
usually performed experiments the electric field can not be directly controlled, we investigate 
the crystal size dependence of the 5-electron LDOS symmetry. It is found that the symmetry 
changes from spherical to torus-like by increasing the nanocrystal radius. 



1. Introduction 

Wave-function mapping in semiconductor quantum dots (QD) has recently attracted much 
interest since it serves as the ultimate tool to study the electronic structure of those dots 
[El S @I- Knowing the actual shape of the electronic densities contributes to a better 
understanding of the QD's electronic structure. This knowledge is crucial with respect to the 
possible importance of semiconductor QDs as the ultimate building blocks of optoelectronic 
and nanoelectronic devices. 

Next to various experimental techniques available for different dot types, recent 
scanning-tunneling-microscopy (STM) measurements also allow an imaging of electronic 
densities in colloidal nanocrystals [jlp. Those nanocrystals are nearly spherical in shape 
resulting in atomic-like symmetries and degeneracies of the electronic states [Q, 
Nevertheless certain experiments [[[[] show densities with a torus-like symmetry. In the 
regime where electrons tunnel through neutral nanocrystals we could attribute this broken 
symmetry previously to the electric field induced by the applied STM voltage. Although 
scanning-tunneling-spectroscopy (STS) experiments on colloidal nanocrystals are possible in 
this transport regime [[7p, in the usually performed experiments the nanocrystal gets charged 
by increasing the STM voltage. Therefore we study in this work the effect of the STM voltage 
on the electronic states in the regime of tunneling through already charged nanocrystals. 
Nevertheless the results of the earlier published single-particle calculation are also of great 
importance for this regime and are therefore reviewed in the first part of this paper. 

We present a calculation of the many-particle conduction band (CB) states within a 
particle-in-a-sphere model taking into account the electric field due to the applied STM 
voltage. We mainly concentrate in this work on a Stark effect induced reordering of states for 
five electrons occupying the nanocrystal. It will be shown that this reordering corresponds to 
a change of the LDOS symmetry which in principal is observable in a wave-function mapping 
experiment. 
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2. Model 



A sketch of the experimental setup in a STS experiment on colloidal InAs nanocrystal 
quantum dots is shown in Fig. [I]. To obtain a tunnel spectrum the tip is positioned above 
a single nanocrystal attached to the substrate via hexane dithiol molecules. Keeping the tip- 
crystal distance constant the differential conductance as a function of applied voltage shows 
sharp peaks [|], |8], ^p. For a detailed understanding of the experimental data it is thus crucial 
to know the discrete energy spectrum of the QD since the obtained peak positions U pea k are 
directly related to the electronic dot structure [@] : 

eU peak (N, N+l) = -y[E{N + 1, a) - E(N, (3)] (1) 

E(N,a = K» = 5> ei + VZ*({m}) (2) 

i 

where the pre-factor 7 depends on the capacitive electrostatic geometry. For a very 
asymmetric tip-dot dot-substrate capacity distribution 7 is close to one The total energy 
E(N,a = {ni}) of excitation level a with N = J^i^i electrons, where rij denotes the 
occupation number of state i, is written as a sum of occupied single particle levels with 
energies and the total charging energy Vq 1 . Hereby Vq 1 includes both, the direct Coulomb 
interaction and all correlations. 




Figure 1. Scanning tunneling spectroscopy of a single InAs nanocrystal. The InAs 
nanocrystal with a typical radius of a few nanometers are linked to a gold substrate by hexane 
dithiol molecules (DT). Trioctylphosphin (TOP) molecules form a ligand shell around the 
nanocrystal. At 4.2K the tunnel current is measured as a function of the applied voltage U 
between tip and substrate. 



In order to obtain the discretised energy levels (g) of the CB electrons we use a single- 
band envelope wave-function approximation. The confinement due to the finite crystal size is 
modeled by a spherical potential well with finite depth [10[]. Depending on the crystal radius 
and the material constants of the used semiconductor those QDs can form a quite strong 
confinement. This confinement can lead to a size quantization of the electronic energies in the 
order of the semiconductor gap. Hence non-parabolicity effects of the CB have to be taken 
into account. Accounting for this effect we use an energy dependent effective mass approach 
pTTp with 

m*{E) =m*(0)[l + E/E g ] (3) 

where m*(0) is the bottom CB effective mass (0.0239 m e in InAs) and E g the bulk energy 
gap. This approach has proven to by quite successful in reproducing the energy gap between 
the single particle ground and first excited state in InAs nanocrystals [Bp. In the STS setup 
of Fig. U with typical voltages up to U ~ 2V \^ |]] applied on a tip-substrate distance of 
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a few nanometers, the QD is exposed to a considerable electric field. Therefore the model 
Hamiltonian for a single CB electron reads 

H = -^J^yV 2 + V e(r - R,) - e$ e (f) (4) 

with $ e (r) being the electrostatic potential. The confinement potential is described by a 
Heaviside function 0(r — Ri) with the potential well depth Vo- 

The electrostatic potential & e (r) is obtained from a realistic modeling of the electrostatic 
environment in the experimental setup of Fig. [I]. While a STM tip has to terminate in a single 
atom in order to achieve atomic resolution, the macroscopic tip size is usually about one 
order in magnitude bigger than the here studied nanocrystals tJT2|]. Other than a macroscopic 
metallic tip, a single terminating atom is not able to substantially focus the electric field. 
Over the nanocrystal size, we can therefore assume the field between tip and substrate to be 
homogeneous in the absence of the QD. Since these nanocrystals are typically surrounded 
by ligands with a quite different relative dielectric constant compared to the semiconducting 
crystal we model the QD as a jacketed dielectric sphere. Extending the textbook calculation 
of a dielectric sphere placed in a homogeneous electric field £h om JT3| ] to such a structure 
leads to a potential inside the QD of 

$(r6) = 9e 2 £ hom rcos9 

H ' ^ 2el + e l£2 + 4e 2 + 2e x + 2(§ ffae, +e 2 -e 1 -e 2 2 ] K) 

with the relative dielectric constants e\ and e 2 of the nanocrystal and the ligand shell, 
respectively. As in the case of a dielectric sphere without a shell JT3| ] the core potential $ e is 
still the potential of a homogeneous field. The field Shom occurring in (|5|) is not directly 
accessible. It is obtained by equating the voltage U with the potential drop between tip 
and substrate of the inhomogeneous field outside the QD. Knowing the electrostatic potential 
inside and also outside the nanocrystal the voltage drop over both tunneling barriers can be 
obtained by plotting the potential along z-direction as shown in Fig. [2| The pre-factor 7 




Figure 2. Electrostatic potential between STM tip and gold substrate divided by the applied 
tip-substrate voltage U plotted along the tip direction. The potential is calculated for an InAs 
nanocrystal of 3.2nm in radius and a tip-crystal distance of lnm. The substrate is on the left 
and the tip starts on the very right. The dielectric constants ei = 15.15 and 62 = 2.1 [|l5| 
have been used for the InAs dot and the ligand shell, respectively. 

of equation (jT)) can also be extracted from the electrostatic potential by simply relating the 
potential drop between tip and QD-center to the drop between QD-center and substrate. 

For nanocrystals charged by iV CB electrons Coulomb interaction between those 
electrons will be present such that the many particle Hamiltonian can be written as 



N 



2m 



-V 2 i -V Q Q{r i -R 1 ) + e$ e 



A' 



+ E 



(6) 
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with e\ being the dielectric constant of the semiconducting nanocrystal. Due to the finite 
potential well depth the electronic wave-functions will to some extend leak out into the ligand 
shell. The much smaller dielectric constant in the ligand shell leads to an enhanced Coulomb 
interaction in this outer region compared to inside the QD. Therefore the Coulomb energy 
will be underestimated by this model Hamiltonian. In fact we found the charging energy for 
InAs nanocrystals about a factor two smaller than the experimentally observed values [|J]. To 
account for this effect one has to replace the Coulomb operator in @ by the proper Green's 
function of a dielectric sphere |TO, [TTQ. Instead we use in this work e% as a fitting parameter to 
match the experimentally found charging energy. Mirror charges induced in the metallic tip 
and substrate are also neglected. 

3. Single-particle calculations 

The single-particle Hamiltonian (Q) without electrostatic potential separates in an angular and 
a radial part where the angular part is solved by spherical harmonics. The radial Schrodinger 
equation is solved by spherical Bessel functions ji inside the well and spherical Hankel 
functions hi outside [|T8[]. The continuity conditions at the potential step lead to a set of 
transcendental equations determining the energy levels 

ahtiipRt) [lji-i{aRi) -(1+ l)j l+1 (aRi)] 

= mi{aRi) [lhi-i(ipRx) - (/ + l)h l+1 (i(3R 1 )] (7) 

with a = \ / 2m*E/h and j3 = \j2m*(V — E)/Ti. For small nanocrystals where the energy 
levels are in the order of the bulk energy gap the energy dependent mass of equation @ can 
be directly inserted into the transcendental equation ([7[). Like in the hydrogen atom the single- 
particle ground-state has a s-type wave-function (hereafter referred to the lS e state, where the 
subscript e denotes an electron rather than a hole state). Other than in hydrogen the allowed 
orbital quantum numbers are not restricted by the principal quantum number. Hence the first 
excited state has the quantum numbers n — 1 and / = 1 which we will call the lP e state. 
Owing to the spherical harmonics, this state is threefold degenerate in the three magnetic 
quantum numbers m=-l, and 1 as shown in Fig. |3]a. 

Knowing the single-particle states without an electric field allows now to calculate the 
Stark effect either by an exact diagonalization or by perturbation theory. In the case of InAs 
nanocrystals the results of a calculation in second order perturbation theory has proven to be 
sufficient Although the lS e state also shows a Stark effect, we concentrate here on the 
first excited lP e state. In contrast to hydrogen, the first excited state does not show a linear 
Stark effect, owing to the lack of s-p degeneracy in such step-like spherical potential. The lP e 
degeneracy is lifted by the quadratic Stark effect, such that the energy of the lP e (m ± 1) wave- 
functions oriented perpendicular to the electric field are lowered compared to the lP e (m = 0) 
wave-function oriented along the field (see Fig. |3]b). 

Using the material parameters of a 3.2 nm InAs crystal leads to a Stark splitting of about 
15 meV (see Fig. ^p) at an applied voltage of 1.4 V which corresponds to the experimental 
situation in [[]]]. This Stark- induced energy splitting is experimentally resolvable and allows 
at an appropriate voltage tunneling into the energetically lower lP e (m ± 1) states without 
tunneling through the lP e (m = 0) state. The qualitatively different density distributions of 
those split states (see Fig. |]d), namely torus-like for the lP e (m ± 1) and spherical for the 
superposition of all three p-type orbitals, are observed in a STM experiment by Millo et al. 

Therefore, the obtained Stark-induced degeneracy lifting of the first excited lP e state 
serves as a possible explanation for the experimentally observed mapping of the lP e (m ± 1) 
wave-functions only. 
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Figure 3. a) Sketch of ground and first excited state without electric field, b) Single particle 
energy levels with an electric field applied along z-direction. The lP e degeneracy is partly 
lifted, c) Stark splitting of the lP e level in a 3.2 nm radius InAs nanocrystal as a function of the 
applied tip-substrate voltage. For this calculation a confining potential well depth of Vo = 3 
eV has been used jl9|]. d) Electronic density of state lP e (m = 0) (top) and lP e (m ± 1) 
(bottom) viewed along the applied electric field. 



4. Many-particle calculation 

In the usually performed experiments the p-type orbitals are available for tunneling only if the 
nanocrystal is already occupied by at least two electrons. Coulomb interaction between the 
CB electrons might be important and therefore the question arises if the single particle results 
of last section concerning the LDOS symmetry are still valid in this transport regime. 

Solutions of the many-particle Schrodinger equation belonging to Hamiltonian @ were 
found by an exact diagonalization procedure in the basis of wave-functions solving the many- 
particle Schrodinger equation without Coulomb interaction and without electrostatic potential. 
In order to keep the obtained matrices as small as possible we used the L z and S z symmetry 
of Hamiltonian (Q) by selecting the needed basis functions. Furthermore the used basis is 
terminated by an energy cutoff, meaning that only Slater determinants with an energy below 
some threshold are used. The accuracy of a calculated energy level is estimated by the relative 
difference to the energy obtained by using a basis of just half the size. For all many-particle 
calculations presented in this work the relative error in energy is smaller than 1.5 • 1CT 3 . 

4.1. Channel with 3 and 4 electrons 

In the case of the first p-type channel a third electron does tunnel through the QD which is 
already occupied by a full l.S'e shell. Owing to the Stark-induced degeneracy lifting the third 
electron will tunnel through the energetically lower lP e {m ± 1) state. Therefore the leading 
term of the 3-electron ground-state corresponds to the configuration shown in Fig. |]a. Due to 
the small crystal size and the high dielectric constant (ei n A s = 15.15) in those semiconducting 
nanocrystals the Coulomb energy is usually small compared to the kinetic energy (about 100 
meV compared to 320 meV in a 3.2 nm InAs dot [Q]). Therefore the LDOS will be dominated 
by this leading configuration and hence torus-like in symmetry. 

Since the lP e (m ± 1) states are still two-fold degenerate the fourth electron will also 
tunnel through one of those energetically lower states. In order to gain exchange energy both 
p-electrons will align their spins (Hund's rule) which is the configuration shown in Fig. 0b. 
The LDOS symmetry of this 4-electron ground-state will again be torus-like. 
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Figure 4. (a) Leading ground-state configuration for three electrons, (b) Ground-state 
configuration according to Hund's rule for 4 electrons occupying the nanocrystal. (c) Two 
possible ground-state configurations for five electrons. 



4.2. Channel with 5 electrons 

A more interesting situation arises for the fifth electron tunneling through the nanocrystal, 
where a competition between Stark energy on one hand and exchange energy on the other 
hand arises. Depending on how strong the splitting of the lP e states is, two different ground- 
states are possible. For a small splitting the configuration with all three p-type orbitals, 
lP e (m = ±1) and lP e (m = 0), occupied each by one electron with their spins aligned will be 
favored (see configuration A in Fig. In this configuration the fifth electron has to pay some 
Stark energy but gains exchange energy. On the other hand if the splitting becomes too big it 
is energetically more favorable for the fifth electron to occupy also a lP e (m± 1) state, thereby 
saving Stark energy. Due to the necessary spin flip however, it has to pay exchange energy 
(see configuration B in Fig. [|c). This competition between Stark and exchange energy leads 
to a ground-state crossing with increasing electric field strength (see Fig. |5|). The interesting 




Figure 5. Ground and first excited state energies of five electrons versus the electric field 
strength inside a 3.2 nm radius InAs nanocrystal. The found ground-state crossing corresponds 
to a change in the LDOS symmetry. For this 3.2 nm dot the 5-electron channel is available 
at a tip-substrate voltage of about 1.5 V leading to an electric field strength indicated by the 
vertical arrow. 

point is that this crossing also corresponds to a change in the LDOS symmetry from spherical 
to torus -like by increasing the electric field. In contrast to the 3- and 4-electron channel the 
LDOS symmetry of the 5-electron channel can be controlled by the electric field applied to 
the nanocrystal. 

Unfortunately the electric field strength applied to the nanocrystal is determined by the 
corresponding peak position found in a STS experiment and, therefore, it is experimentally 
not straight forward to switch the LDOS symmetry forth and back between torus-like and 
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spherical in a wave-function mapping experiment. For a 3.2 nm InAs nanocrystal the 5- 
electron channel is available at a tip-substrate voltage of about 1.5 V leading to an electric field 
of about 0.1 V/nm. As indicated by the vertical arrow in Fig. ^| the LDOS symmetry of the 
ground-state is at this field strength still spherical. On the other hand the STS-peak positions 
depend on the QD radius such that we now stress the question how the LDOS symmetry of 
the 5-electron ground-state changes with the nanocrystal radius. 

In order to answer the question how the ground-state symmetry depends on the crystal 
size the scaling behavior of the Coulomb operator versus the electrostatic potential in (Q) with 
respect to the dot radius is studied. Whereas it is clear that the Coulomb operator scales with 
i?f 1 the scaling of the electrostatic potential is not easily foreseen. As shown in equation (0) 
the electrostatic potential $ e oc £hom(Ri) • R\ scales linearly with the dot radius and electric 
field strength. This field also depends on the dot size, but other than in a plain capacitor it 
scales roughly with Shom oc U(R\) ■ i?f ' 4 . The reason for this scaling behavior is mainly 
the fact that the tip-crystal distance is kept constant while scaling the crystal size. Last but 
not least the applied tip-substrate voltage depends on the energy needed to add a further 
electron to the crystal which is again a function of the dot radius. In an infinite potential 
well the single-particle energy levels scale with i?f 2 but due to the finiteness of the studied 
potential well and the fact that the effective mass increases with increasing energy, leads to a 
scaling of roughly U oc iff . Putting all together we find that the electrostatic potential scales 
with $ e oc i?f °' 4 . Therefore the LDOS symmetry changes from spherical to torus-like with 
increasing the crystal radius, since Stark energy becomes in bigger crystals more important 
than exchange energy. 

Especially the scaling of the electrostatic potential is not straight forward and therefore it 
is necessary to check this result by a full calculation. To this end, the charging energy needed 
to add the fifth electron to the QD has to be calculated in a first step. As shown in equation 
(0) this charging energy is the energy difference between the 5- and 4-electron ground-states. 
In a second step the tip-substrate voltage needed to open this 5-electron channel has to be 
calculated. As shown in equation (jl[) this voltage is found by multiplying the charging energy 
with the pre-factor 7, obtained from the electrostatic potential drop along z-direction (see 
Section 0). Knowing the applied tip-substrate voltage, the electric field strength in the QD 
can in a last step be calculated by equation (Q). Since the charging energy calculated in the 
first step also depends on the electric field, the whole cycle is repeated until self-consistency 
is obtained. Now knowing the electric field strength we can go back into Fig. |5| and determine 
the ground- state configuration for the considered QD size and therefore determine the LDOS 
symmetry. We have done this calculation for six InAs crystals with radii between about 2 and 
6 nm and plotted the energy difference between configuration A and B versus the dot radius 
in Fig. In this plot a number smaller than zero corresponds to a spherical and a number 
bigger than zero to a torus-like ground-state symmetry. As already predicted by the scaling 
considerations we found that the LDOS symmetry changes from spherical to torus-like by 
increasing the crystal radius. Using the material constants for InAs we found this transition to 
happen at a dot radius of about 4 nm. 

5. Conclusion 

We calculated the Stark effect on many-particle wave-functions in InAs nanocrystals charged 
with up to five electrons. Stark effect and Coulomb interactions have been fully included 
by an exact diagonalization procedure. We found that the LDOS, due to the Stark-induced 
degeneracy lifting of the lP e states, of the 3- and 4-electron ground-states are torus-like in 
shape. For 5 electrons however, a competition between Stark energy and exchange energy 
leads to a ground-state crossing with increasing field. Since the electric field strength can 
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Figure 6. The difference in energy between the 5-electron ground-state and first excited state 
is shown as a function of nanocrystal radius. For an energy difference greater than zero, 
configuration B shown in Fig. Qc is the new ground-state which has a torus-like LDOS 
symmetry. 

not be directly controlled in the usual experiments, we studied the crystal size dependence of 
the 5-electron ground-state symmetry and found a transition from spherical to torus-like by 
increasing the dot radius. This transition should be observable in a wave-function mapping 
experiment. 
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